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ABSTRACT 


The  lethal  volume  of  a  continuous  rod  warhead  is  defined  as 
that  volume  in  which  a  target  reference  point  must  be  located  at 
burst  if  the  target  is  to  be  subjected  to  a  lethal  cut  from  the 
rod.  The  lethal  volume  is  a  function  of  the  characteristics  of 
the  target,  warhead,  and  intercept  dynamics. 

A  continuous  rod  warhead  subjects  a  target  to  lethal  damage 
by  ejecting  a  Icop  of  metal  laterally  to  cut  a  structural  member 
in  the  target.  As  the  loop  is  ejected  it  expands  in  a  circle 
until  the  radius  of  breakage  is  reached. 

The  target  is  assumed  to  be  a  single,  linear,  rigid,  and 
non-rotating  structural  member.  The  lethal  volume  is  obtained 
for  this  target  by  examining  the  relative  motion  of  the  target 
and  of  the  rod  with  respect  to  the  missile  carrying  the  warhead. 
Under  the  assumption  of  constant  rod  velocity,  the  lethal  volume 
is  shown  to  be  composed  of  linear  segments  joining  two  identical 
right  elliptic  cones. 

Expressions  for  the  size  and  probsblllty  content  of  the  lethal 
volume  are  then  obtained  for  the  constant  rod  velocity  case.  A 
numerical  schema  is  discussed  with  which  the  probability  content 
of  the  lethal  volume  can  be  evaluated. 

Parameter  variation  studies,  which  can  be  performed  with  the 
lethal  volume  nodel,  are  next  discussed.  As  examples,  the  varia¬ 
tion  in  the  size  and  probability  content  of  the  lethal  volume  is 
computed  for  several  combinations  of  the  input  parameters. 

The  thesis  is  closed  with  two  proposed  schemes  for  improving 
the  model.  First,  a  means  of  improving  the  target  model  such  thst 


it  bears  a  closer  resemblance  to  an  airplane  Is  discussed.  Fina 
ly,  a  means  of  accounting  for  variable  rod  velocity  is  given. 


CHAPTER  I 


INTRODUCTION 

A.  LETHAL  VOLUMES 

The  lethal  volume  of  a  warhead  Is  that  volume  in  which  a 
target  reference  point  must  be  located  at  the  time  of  burst  If  it 
is  to  be  subjected  to  some  lethal  effect  from  the  warhead.  Thus  a 
necessary  and  sufficient  condition  for  target  kill  is  that  the 
reference  point  be  located  in  the  lethal  volume  at  the  time  of  burst. 
The  location  of  the  lethal  volume  is  usually  considered  relative  to 
the  center  of  the  warhcd,  but  is  occasionally  considered  relative 
to  some  other  point  such  as  a  point  on  the  target.  Here,  the  former 
procedure  will  be  followed. 

As  an  elementary  example  of  a  lethal  volume,  consider  a  warhead 
which  destroys  all  point  targets  within  a  distance  R  of  the  center 
of  the  warhead  at  the  time  of  burst.  The  lethal  volume  would  then 
be  a  sphere  of  radius  R  centered  about  the  warhead.  The  distance 
R  is  called  the  lethal  radius  of  the  warhead.  The  lethal  radius  is 
primarily  a  function  of  the  warhead  and  target  characteristics  but 
is  also  influenced  by  environmental  factors  such  as  atmospheric 
density. 

The  usefulness  of  the  concept  of  a  lethal  volume  lies  largely 
in  its  relation  to  the  kill  probability  of  an  operable  missile.  A 
missile  is  said  to  be  operable  unless  it  fails  because  of  weapon 
system  malfunctions.  The  kill  probability  of  an  operable  missile, 

P,  is  defined  as  the  probability  that,  under  a  given  set  of  target, 
missile,  warhead  and  intercept  conditions,  an  operable  missile  will 
result  in  a  target  kill. 
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If  f(x,y,z,)  is  the  probability  density  function  describing 
the  location  of  the  target  reference  point  relative  to  the  warhead 
at  the  time  of  burst,  then  the  kill  probability  of  an  operable 
missile  is  given  by, 


p  ■  J  J  J  f(x,y,z)dxdydz  , 


where  V  is  the  lethal  volume.  If, 

„  *  \ _ 1  .  Cx2  +y2 +z3  ] 


f(x,y,z) 


(2nj)z 


describes  the  distribution  of  point  targets  about  the  center  of 
the  warhead  for  the  spherical  lethal  volume  considered  earlier, 
then  a ,  a  a  q  a 

P . r r r _i  t* *r « . 

J  o3  ,  . 

If  P(W  <  R)denotes  the  probability  that  v<H,then  P  =  pf X  — 2  < 

O 

This  is  simply  P  (x2  <  R2 /a2') 

>  8  “  ' 


where  ^  has  a  chi-square  distribution  with  three  degrees  of 
freedom. *  This  probability  can  be  found  from  a  table  of  the 
chi-square  cumulative  distribution. 

In  practice,  it  is  rarely  possible  to  evaluate  P  in  terms  of 
tabulated  functions.  For  this  reason,  numerical  integration  of  the 
integral  expression  for  P  must  normally  be  used. 

Once  the  relationship  between  the  lethal  volume  and  the  kill 
probability  has  been  established,  system  parameter  tradeoffs  can 
be  made.  To  perform  tradeoffs,  system  parameters  are  relsted  to 
the  lethal  radius  of  the  warhead  and  the  guidance  accuracy  of  the 


^H.  W.  Lilliefors,  *'A  Hand -Computation  Determination  of  Kill 
Probability  for  Weapons  Having  Spherical  Lethal  Volume",  OPERATIONS 
RESEARCH  JOURNAL,  Vol.  5,  No.  3,  June  1957,  pp ,  416-421. 
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missile.  These  parameters  can  be  varied  and  their  effect  on  kill 
probability  observed.  If  the  model  is  sensitive  to  environmental 
factors,  their  effects  on  kill  probability  can  also  be  investigated. 
On  the  other  side,  studying  the  effects  of  hardening  the  target  or 
forcing  the  burst  to  occur  closer  to  the  target  is  a  problem  fre¬ 
quently  investigated.  In  all  cases  significant  changes  in  the  ki'.  1 
probability  are  sought  for  a  minimum  of  cost. 

In  addition  to  the  kill  probability,  another  measure  of  warhead 
effectiveness  which  is  often  used  is  the  size  of  the  lethal  volume. 
Since  the  kill  probability  is  obtained  by  integrating  over  the 
lethal  volume,  the  variation  of  the  size  of  the  lethal  volume  can 
be  used  as  a  crude  measure  of  the  variation  in  kill  probability. 

ihe  size  of  the  lethal  volume  as  a  measure  of  effectiveness 
is  useful  since  it  is  easier  and  cheaper  to  compute  than  the 
probability  of  kill.  For  instance,  the  IBM  7090  computer  program 
for  the  lethal  volume  of  a  continuous  rod  warhead  based  on  this 

thesis  computes  the  size  of  the  volume  at  less  than  one -fourth 

2 

the  cost  of  the  calculation  of  kill  probability.  For  this  reason, 
the  size  of  the  lethal  volume  is  preferred  as  a  measure  of  warhead 
effectiveness  when  the  size  of  the  volume  adequately  reflects  the 
parameters  under  study. 

B.  CONTINUOUS  ROD  WARHEADS 

To  form  a  continuous  rod  warhead,  a  single  long  continuous 
metal  rod  is  folded  onto  the  warhead  over  a  high  explosive  charge. 

'  2  "  " 

Vitro  Laboratories  Program  R-16,  Silver  Spring,  Maryland, 


1964. 
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Upon  detonation  the  rod  expands  in  a  circle  normal  to  and  centered 
at  the  missile  axis,  as  shown  in  figure  I  A.  The  circle  of  metal 
expands  as  a  continuous  rod  until  the  circumference  of  the  circle 
is  approximately  equal  to  the  length  of  the  unfolded  rod.  Then 
the  rod  breaks  randomly  into  discrete  lengths.  The  radius  of  the 

i 

1 

circle  of  metal  at  the  time  of  breakup  is  called  the  maximum  opening 
radius  of  the  rod  and  will  be  denoted  by  R^. 

The  speed,  thickness,  and  mass  of  the  rod  are  such  that  if  the 
rod  strikes  a  structural  member  of  an  aircraft  the  member  will  be 
severed.  Such  damage  will  usually  cause  flight  failure  unless  the 
member  is  cut  on  an  extremity  such  as  a  wing  tip.  After  the  rod 
reaches  its  maximum  opening  radius,  the  effectiveness  of  the  rod 
falls  off  rapidly  and  can  be  ignored  in  kill  probability  estimates. 
(See  figure  IS) . 

The  lethal  volume  of  a  continuous  rod  warhead  is  then  defined 

as  that  volume  in  which  the  target  reference  point  must  be  located 

at  burst  if  the  target  is  to  be  lethally  cut  by  the  rod  before  the 

3 

rod  reaches  Its  maximum  opening  radius  R^. 

C.  PROBLEM  STATEMENT 

The  major  problem  addressed  in  this  thesis  was  to  mathematically 

a 

describe  the  set  of  locations  of  a  target  reference  point  relative 
to  the  warhead  at  the  time  of  warhead  detonation,  such  that  the 
subsequent  motion  of  the  target  and  the  rod  would  resulL  in  the 
rod  cutting  the  target.  This  set  was  defined  as  the  lethal  volume 
of  a  continuous  rod  warhead. 

^The  contents  of  the  section  were  paraphrased  from  M.C.Waddcl, 
Surface-To-Air  Guided  Missile  Systems,  Applied  Physics  Laboratory 
TG  396,  March  1961,  pp.  37-38. 


MISSILE 

VELOCITY 


Figure  IA  ROD  EXPANSION  RELATIVE  TO  THE  GROUND 
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Since  a  necessary  and  sufficient  condition  for  the  target  to 
be  cut  by  the  rod  (target  kill)  is  that  the  reference  point  be  in 
the  lethal  volume  at  detonation,  the  probability  of  kill  is  the 
probability  that  the  reference  point  is  in  the  lethal  volume  at 
the  time  of  warhead  burst.  A  subsequent  problem,  then,  was  to 
obtain  the  probability  content  of  the  lethal  volume  and  the 
variation  in  this  probability  content  with  model  inputs. 

D.  IMPORTANCE  OF  STUDY 

To  date  there  has  been  very  little  effort  devoted  to  an 
analytical  description  of  the  lethal  volume  of  a  continuous  rod 
warhead.  Nearly  all  effort  in  this  area  has  been  devoted  to 
large  scale  simulations  in  which  specific  targets,  missiles,  and 
warheads  are  "flown"  against  each  other  under  a  given  set  of 
intercept  or  end  game  conditions.  Using  Monte  Carlo  techniques, 
a  specific  end  game  is  run  repeatedly  until  a  probability  of  kill 
is  obtained.  The  simulation  technique  has  two  major  disadvantages. 
Firat,  it  is  inherently  expensive.  Secondly,  it  fails  to  give  the 
analyst  an  adequate  "feel"  for  the  relative  importance  of  the 
variables  under  study.  Because  of  the  tremendous  costs  associated 
with  the  simulation  approach,  tradeoff  studies  among  the  inputs 
oust  be  limited  in  scope. 

On  the  other  end  of  the  spectrum  of  the  tools  that  exist  for 
the  analysis  of  continuous  rod  warheads,  are  some  elementary 
analytical  models.  Those  models  are  accurate  only  for  the  simplest 
intercept  geometries  such  as  those  in  which  the  velocity  vectors 
of  both  the  target  and  the  missile  are  along  a  straight  line.** 


4 Ibid .  pp,  51-53. 
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The  major  disadvantage  of  these  models  is  their  insensitivity  to 
many  of  the  parameters  which  are  known  to  significantly  affect 
the  kill  probability  associated  with  a  given  end  game. 

The  lethal  volume  model  developed  here  was  designed  to 
bridge  the  gap  between  the  simulation  technique  and  simplified 
model;  i.e.,  the  lethal  volume  was  intended  to  be  sensitive  to 
many  of  the  variables  which  the  simplified  model  is  not  and  at 
the  same  time  be  substantially  cheaper  to  use  than  the  simulation 
technique.  Because  of  the  large  decrease  in  costs,  tradeoff 
studies  performed  with  the  lethal  volume  model  can  be  much  more 
complete  than  is  currently  possible  with  the  simulation  technique. 
For  this  reason,  the  lethal  volume  model  is  believed  to  be  a 
significant  new  tool  for  the  analysis  of  continuous  rod  warhead 
problems . 

E.  ORGANIZATION  OF  THESIS 

Chapter  II  of  the  thesis  contains  the  derivation  of  the 
lethal  volume  of  a  continuous  rod  warhead  against  a  single,  linear, 
non-rotating,  and  rigid  line  segment.  The  line  segment  is  assumed 
to  represent  the  fuselage  or  wing  of  an  aircraft.  The  lethal 
volume  is  obtained  for  this  target  by  examining  the  relative 
motion  of  the  target  and  the  rod  with  respect  to  the  missile 
carrying  the  warhead.  Chapter  II  is  closed  with  examples  of  the 
shape  of  the  lethal  volume  as  a  function  of  the  relative  location 
of  the  target  and  the  missile  carrying  the  continuous  rod  warhead 
at  burst. 

Chspter  III  contains  expressions  for  the  size  and  probability 
content  of  the  lethal  volume.  A  numerical  scheme  for  the  evnlua- 
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tion  of  the  expression  of  the  probability  content  of  the  lethal 
volume  is  discussed  and  illustrated. 

Chapter  IV  is  devoted  to  an  examination  of  the  uses  of  the 
expressions  for  the  size  and  probability  content  of  the  lethal 
volume  derived  in  Chapter  III.  These  uses  arc  concerned  with 
the  variation  in  the  size  and  probability  content  with  the 
characteristics  of  the  warhead,  the  missile  carrying  the  warhead, 
and  the  target.  A  sample  study  concludes  Chapter  IV. 

Chapter  V  is  devoted  to  means  of  improving  the  warhead  and 
target  models  used  in  the  derivation.  Particular  attention  is 
given  to  improving  the  target  model  so  that  is  bears  a  closer 
resemblance  to  an  aircraft. 


CHAPTER  II 


DERIVATION 

A.  COORDINATE  SYSTEMS  AND  DEFINITIONS 

In  this  section  two  coordinate  systems  will  be  defined.  One 
system  will  travel  with  the  missile;  therefore,  the  rod  will  remain 
in  a  plane  in  this  coordinate  system.  Because  of  this  fact  the 
equations  for  the  lethal  volume  are  greatly  simplified  in  this 
coordinate  system.  A  second  coordinate  system  will  be  defined 
such  that  it  has  the  same  orientation  as  the  first  but  will  be 
stationary.  The  second  system  will  then  be  used  to  establish  the 
motion  between  a  rigid  structural  member  (target)  and  the  missile. 

Let  M  be  the  point  at  the  center  of  the  warhead.  A  right- 
handed  orthogonal  coordinate  system  with  its  origin  at  M  will  now 
be  established.  First  the  x-ax:ls  is  chosen  such  that  it  lies  along 
the  axis  of  the  missile.  The  z-sxis  is  then  chosen  such  that  it 
is  perpendicular  to  the  x-axis  and  such  that  the  x-z  axes  lie  in 
the  vertical  plane,  with  positive  z  "above”  the  origin.  The  y-axis 
is  chosen  to  complete  a  right-handed  system,  (see  Figure  2).  Let 
I,  J,  and  It  be  unit  vectors  directed  along  the  positive  x,  y,  and 
z  axes  respectively. 

Since  this  coordinate  system  moves  with  the  missile,  the  rod 
will  move  in  the  y-z  plane.  Furthermore,  since  the  rod  expands 
in  a  circle,  the  path  of  the  rod  is  described  by, 

K  -  (R  sin  0)3  ♦  (R  cos  9)E 

where , 

R  is  the  radius  of  the  circle, 

9  is  the  angle  between  R  and  E. 


figure  2  COORDINATE  SYSTEMS 
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Note  that  R  <  RL>  the  maximum  opening  radius. 

Let  G  be  the  origin  of  a  second  coordinate  system,  having 
the  same  orientation  as  the  above  system  but  fixed  relative  to 
the  ground.  This  coordinate  system  will  be  used  only  to  establish 
the  relative  motion  between  M  and  the  target. 

Now  consider  a  rigid  linear  structural  member  representing  a 
major  component  of  an  aircraft  such  as  the  fuselage  or  a  wing. 
Define  B  and  C  as  the  end  points  of  the  largest  segment  on  the 
member  such  that  a  cut  of  the  segment  will  be  lethal  to  the  member. 

Define  8£  as  the  vector,  having  length  equal  to  the  distance 
from  B  to  C  and  having  the  direction  from  B  to  C.  Let  the 
direction  cosines  of  BC  be  1  ^  .  and  L  .  |!}Cl  will  be  referred 

13  3 

y 

to  as  the  effective  length  of  the  member. 

Let  A  be  a  point  on  the  member  between  B  and  C.  The  velocity 
of  A  relative  to  H  will  now  be  established.  Referring  to  figure  2, 
the  vector  MA  is  seen  to  be  given  by: 

MA  ■  CX  ~  CM. 

Differentiating  both  sides  with  respect  to  time,  t»  gives 


d  MA 


d  MA  -  d  GA  -  d  GM 
dt  dt  dt 


Defining  as  the  velocity  of  A  relative  to  M,  denoted 

*k/H  ’  glVC3 

TVh  ’ 

Under  the  assumption  that  there  is  no  rotation  of  the  member 
during  the  period  of  interest,  the  velocity  of  all  points  of  the 
member  will  be  equal  to  .  This  assumption  is  adequate 

because  of  the  extremely  short  time  between  warhead  burst  and  the 


time  the  rod  reaches  its  maximum  opening  radius  R^.  Denoting  this 
common  velocity  by  gives 

rK'V/ 'O-Wo  • 

Let  the  direction  cosines  of  V.  ».  be  m  ,  m  ,  and  m  and  the 


absolute  value  of  be  I 


In  addition. let  the  direction 


cosines  of  Vu/„  be  n  ,  n  ,  and  n  and  let  the  absolute  value  of 


^m/o  be  *  Vg  *  *  11160 

VM/G  "  ^  VM/G  '  ni  1  +  *  VK/&  *  ns  ^  +  '  VM/G  1  ns  E 
and 


va/o  "  1  \/a  1  mi  i^Vo|m»  I4,Vo' 

Substituting  these  values  into  equation  (1) 

-  wOMWv 

4(lVbln.  _l  V&K)* 


In)  j 


where. 


V  T+V  J  *  V  lc 
1  r  » 


v,  Volni 


",  1  k/a '  ~ '  Vo  '  n, 

\  • 1  Vo 1  ‘  1  Vo 1  n« 


B.  PARAMETRIC  EQUATION  OF  THE  LETHAL  VOLUME 

All  positions,  relative  to  M,  at  which  the  point  B  can  be 


located  at  burst  and  the  member  still  be  cut  by  the  rod  between 
B  and  C  will  now  be  determined.  The  location  of  the  point  B  will 
be  defined  such  that  any  point  of  the  member  will  be  cut  by  the  rod 


at  a  radius  R  and  angle  9  .  The  selection  of  the  point  B  as  the 
target  reference  point  is  only  for  convenience  since  any  target 
reference  point  would  suffice. 

The  time  for  the  rod  to  travel  a  distance  R  is  given  by 

t  -  R/Vr 

where, 

VR  is  the  average  velocity  of  rod  to  the  distance  R.  If  A  is 
an  arbitrary  point  on  the  member,  the  distance,  with  respect  to  M, 
that  A  travels  during  this  time  is  I  3^  I  where 


Referring  to  figure  3,  R  is  seen  to  be  given  by 

K  -  13(A)  +  +  of 

where,  B(A)  is  the  set  of  coordinates  of  the  reference  point  B  (at 
the  time  of  burst)  such  that  the  subsequent  motion  of  the  member 

and  the  rod  will  result  in  the  rod  being  cut  at  the  point  A. 

—  —  —  —  BT 

Since  BA  lies  along  BC,  BA  «XBC,  where  X  ■ - 

BC 

and  0<X  <  1.  Substituting  and  solving  for  B(A)  gives 

B(A)  -  R  -  dj[  -  XBC  (5) 


when  X  ■  0, 

5(B) 

and  when  X  ■ 

5(C)  - 

s  -  r;  -  is  . 

Therefore, 

B(C)  - 

B(B)  -  BC 

(6) 

5(A)  - 

B(B)  -  XBC 

and 

B(A)  - 

B(B)  -  X[B(B)  -  B(C)] 

(7) 

Figure  3.  PATH  OF  TARGET  TO  BE  CUT  AT 
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Two  important  observations  concerning  the  volume  generated  by 
B(A)  can  be  made  from  equations  (?)  and  (7) .  Equation  (7)  shows 
that  the  volume  generated  by  B(A?  is  composed  of  linear  segments 
joining  the  surfaces  generated  by  B(B)  and  B(C) .  Equation  (6)  shows 
that  the  surface  generated  by  B(C)  is  identical  to  that  generated 
by  B(B)  but  is  displaced  by  the  vector  -BC.  Hence,  describing  the 
lethal  volume  generated  by  B(A)  is  essentially  reduced  to  describing 
the  surface  generated  by  B(B) ;  i.e.,  the  surface  associated  with  cuts 


of  the  point  B. 

Substituting  for  R,d7  ,  and  BC  in  equation  (7)  gives, 

A 

1(A)  •  [(5-)  V  .  A|55  U  J:  *  [a  cos  9  4  (|-)  -  A  I  BC  | 

♦  [r  sin  8  ♦  (|-)  V  -  A'  13  |  iji  -  [v  yA,  .  (S) 

Equation  (8)  is  the  parametric  equation  of  the  lethal  volume. 


When  evaluated  for  all  R,  0,  andXsuch  that  0  <  R  <  R^,  0  <  9  ^  2n, 
and  0  <  X  <  1  every  point  at  which  the  reference  B  can  be 


located  at  burst  and  the  target  be  cut  between  B  and  C  is  generated. 
The  volume  generated  by  B(A)  will  be  denoted  V. 

D.  ELIMINATION  OF  PARAMETERS 


Although  equation  (8)  can  be  used  to  generate  the  lechal 

volume  V,  the  elimination  of  the  parameters  R  end  0  will  lead  to 

a  more  efficient  computation.  Since  no  equation  exists  for  reliably 

expressing  Vjj  as  a  function  of  R,  a  completely  closed  form  solution 

is  impossible  to  obtain.  However,  a  closed  form  solution  can  be 

obtained  for  any  given  R.  In  addition,  under  the  assumption  that 

V  is  constant,  a  closed  form  solution  for  the  entire  volume  will 
K 
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be  obtained. 

As  pointed  out  in  the  preceding  section,  the  problem  of 
describing  V  is  essentially  equivalent  to  that  of  describing 


the  surface  generated  by  B(B) .  To  find  B(B) ,  X  =  0  is 

substituted 

in  equation  (8)  giving, 

■B(s)  -  [xB,  yB>  *B]  *  (jjj)  v.  T  *  [R  cos  e  4  Cv^.'  vr  j3 

♦  [r  sin  9  +  (J-)  y 

or 

*B  ’  Q~)  T, 

(9A) 

yg  -  R  cos  9  +  (^~) 

(9B) 

*3  *  R  sin  0  +  (^)  Vt 

(9C) 

Eliminating  9  from  equations (9B)  and  (9C)  gives, 
o  o  2V  R8  r, 

*B4  rO  • 

R  VR 


But  from  equation  (9A)  is  obtained, 

,  if  V  /  0 


R  ■  v 

«  vb  ^ 


Substituting  gives, 

V3  ♦  v3  -  vl 


2V  2V 

+  *b  -  nf  ¥b  -  nr  XBSB  *  0 


(10A) 


where, 

H  •  (V*  ♦  V3  -  v?)/v*  (IOC) 

sr  *  n  * 

B(B)  -  [IB,  yB,  IB] 

Cb(b)]t  is  the  transpose  of  3(B) 

If  M  denotes  the  matrix  of  coefficients  in  equation  (10B), 

then 

S(B)  -  [B(B)]  M  [B(B)]T  (11) 

Where  S(B)  =  0  is  the  equation  of  the  surface  generated  by  B(B); 
i.e.,  the  equation  whose  locus  is  the  surface  associated  with  cuts 
of  the  point  B. 

C.  REDUCTION  TO  NORMAL  FORM 

In  this  section  the  normal  form  of  S(B)  will  be  obtained. 

First,  the  special  case  when  both  V  and  V^are  zero,  will  be  treated. 
Then  the  canonical  form  for  the  case  when  V  and  V  are  not  both  zero 

•  T 

will  be  obtained. 

If  an  orthogonal  matrix  Q  can  be  found  such  that, 

QHQT  -  D  (12) 

(12A) 


where  D  is  a  diagonal  matrix,  then  the  substitution 

B(B)  -  B'(B)  Q 
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can  be  made  in  equation  (11) ,  giving, 

S(B)  -  b'(B)QM(b'(E)Q)T«  B'(B)y>^T[B'(B)]T  -  3  '( B)D[  B  '(  B)  ]T 

This  will  then  be  the  required  normal  form  of  S(B)  in  terms  of 

E'(B)  -  [Xg,  y',  Zg  ].  The  relationship  between  B(B)  and  can 

be  obtained  from  equation  (12A). 

I 

First,  consider  the  special  case  when  both  V  and  V  are  zero. 

r  * 


Then 


M  - 


5 

V* 

s 


0 

0 


1 

0 


0 

1 


and. 


S(B)  -  +  - 


v2 


2  2 

From  equation  (9C)  the  maximum  value  of  z0  is  seen  to  be  R  . 

Xj 


Dividing  through  by  R  gives 


2  a 
"  yB 


'b 


B 


K8 


'  €59 


(13) 


as  normalized  equation  of  the  surface.  For  a  given  R,  V_  is  a 

R 

constant  and  the  surface  associated  with  '■its  of  the  point  B  is 
seen  to  be  a  right  circular  cone.  If  VDis  assumed  to  be  constant 
for  all  R,  then  the  cone  for  R  »  R^will  contain  the  cones  for  all 
other  values  of  R.  Under  this  assumption,  equation  (13)  takes  the 
form, 


*  14  <$) 
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Next  consider  the  case  when  neither  V  nor  Y  is  fcero.  Following 

7  « 

the  characteristic  root  -  characteristic  vector  scheme  for  deriving 
the  characteristic  roots  were  found  to  be, 


(H  +  l)3  +  h 


H  +  1  + 


T 

_R 

V3 

—  >  0 


H 


2  ^ 

,  (H  +  l)2  «■  h  -j 


(l$A) 


(15B) 


and 


p3  -i 


(15c) 


The  characteristic  vectors  were  found  to  be 


"  t1'  IT 


V  1  V  1 

ITT  V  '  I =p“  V~J  » 

i  >  it 


(15D) 


-  r,  _^T-i 

%  “  L1'  1-d  v  *  l-p  V~J 

St 


a  t 


(l^E) 


and 

Therefore 


%-[o,W]  . 


VI"!  | 

V I9*  t 

LV '  %  >J 

Substituting  Q  in  equation  (12)  gives 


D  - 


Hence  S(B)  has  form 


POO 

1 


P  0 

s 

0  1 


S(B)  -  *  ex'"  -  da  y £ 


B 


B 


where,  c  »  Jp~  and  d  *  -J-o 
i  a 


(15F) 
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The  trana forma  Lion  equations  are  given  bv 


b'(b)  «  b(b)qt 


*B  "  [XB  +  V*"  Cup  )  yB  +  (tfO  Cup  )  Zb]  /  '  qi  ' 


(16A ) 


yB  “  [XB  +  T’  G“5“) 


-J  yB +  Vtv1  Virr)  ZB]  / 1  %  1 

a  i  a 


;B  “  [yB  -  t  2b1  /  ’  %  I 


(163) 


(l6C) 


Substituting  from  equations  (9B)  and  (9C)  into  (16C)  gives 


'  r  1 

ig  »  (R  cos  9  -  R  sin  9^ 


V  +  V* 


for  a  given  R,  the  maximum  value  of  Zg  is  now  desired  in  order  to 

establish  the  general  form  of  S(B)  .  Differentiating  with  respect 

to  9  and  equating  to  zero  gives, 

-7 

tan  9  *  i/ 

Substituting  into  equation  (17)  yields 

R  V  R  R  V 


(*B>  *  — 7 - V"  +  ~5 - • 

B  Max  (V  +  tf)  (7*  +  V*) 


a  T 


*  y 


Therefore  the  normalized  form  of  the  surface  associated  with  cuts 


of  the  point  B  is 


A  A 


*B  ^  yB 


R8  +  ,R^ 


(|)  (§) 


Note  that  when  either  7  or  7  but  not  both  is  zero  the  vector 

y  * 


must  be  modified.  When  7  is  zero,  then 

y 


*B  “  yB 


and  when  7^  »  0, 


then  z 


B  “  XB 


(19A) 

(19B) 


21 


In  both  cases  (z*  )>lax  again  equals  R.  However,  equation  (18)  still 

D 

holds  since  the  formulas  for  c  and  d  remain  unchanged.  If  both  Vy 
and  Vz  are  zero  the  characteristic  values  change  and  equation  (13) 
must  be  used,  i'gain,  the  surface  associated  with  cuts  c;  the 
point  B  is  seen  to  be  a  right  cone  or  more  specifically  in  this 
case  a  right  elliptic  cone.  As  before,  if  VR  is  assumed  to  be 
constant,  then  the  cone  for  R  =  will  contain  all  other  cones 
and  equation  (18)  becomes  a  z  3 

n_ 

v  (Vf(£f 

A  definitive  statement  about  the  lethal  volume  V  can  now  be 

made.  Earlier  the  lethal  volume  was  shown  to  be  composed  of  all 

linear  segments  of  length  |  BC|  joining  two  identical  surfaces 

|  _  /■ 

generated  by  B(B)  and  B(C)  where  I  BC I  was  the  effective  length  of 

the  member.  These  facts  in  conjunction  with  the  foregoing  state¬ 
ments  about  the  surface  generated  by  B(B),  gives  the  following 
theorem: 

Under  the  assumption  of  a  constant  non-zero  rod 
velocity  (V_/0)  and  a  non-zero  closing  velocity 

projected  on  the  missile  axis  (V  /  0) ,  the  lethal 

s 

volume  of  a  continuous  rod  warhead  directed  against 
a  linear,  rigid,  and  non-rotating  structural  member, 
is  composed  of  linear  segments  connecting  the  surfaces 
of  two  right  elliptic  cones.  Furthermore,  these  cones 
are  identical  in  size  and  orientation.  They  are  displaced 
from  each  other  by  a  vector  whose  length  is  equal  to 
the  effective  length  of  the  member  and  whose  direction  is 


given  by  the  orientation  of  the  member  with  respect  to 
the  axis  of  the  missile  carrying  the  warhead. 
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EXAMPLES 

Before  proceeding  further  with  the  derivation  three  simple 
examples  will  be  given.  The  conditions  associated  with  the  examples 
will  be  referred  to  as  end  game  conditions.  End  Game  I  is  the 
simple  case  when  the  velocity  vector  of  the  target  and  the  target 
velocity  vector  of  the  missile  both  lie  along  the  same  line.  In 
End  Game  II  the  member  is  inclined  at  45“  to  the  axis  of  the 
missile.  In  End  Game  III  the  member  is  inclined  at  90®  to  the 
axis  of  the  missile.  Table  I  defines  the  three  end  games  in  detail. 
The  examples  were  treated  parametrically  except  when  complexity 
dictated  otherwise. 


TABLE  I.  Definition  of  End  Games 


Model  Input** 

End  Game  I 

End  Game  II 

End  Game  III 

if 

Rr 

\ 

*L 

\ 

VR  (constant)  iiOO  fps 

iiOO  fps 

\>  V 

l  ,  l  ,  t 
i'  a*  a 

— <V^ 

T  1  °* 

0^  0| 

1  BO  | 

|  1c  1 

lie  | 

Ibc  I 

n  .  m  .  n 
l '  a'  a 

-1,  0,  0 

-fir t  n  — fir 
— 

0,  0,  -1 

n  .  n  .  n 
l*  a'  a 

1,  0,  0 

1,  0,  0 

If  o*  o 

'Vo' 

300  fps 

300  fps 

200  fps 

200  fps 

^These  terms 

are  defined  on 

pp. 11-13 
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In  End  Game  I  both  Vy  and  V2  are  zero  and  equation  (I/*)  holds, 

consequently,  the  surface  associated  with  cuts  of  the  point  B  is  a 

1  V  IfL 

right  circular  cone  of  height  -i - — 


Figure  4  shows  the  lethal  volume  for  this  end  game. 

i 

In  End  Game  II  the  member  is  inclined  at  45°  to  the  axis  of 
the  missile.  The  velocity  vector  of  the  member  and  the  velocity 
vector  of  the  missile  are  aligned  with  the  member  and  the  axis  of 
the  missile,  respectively.  From  equations  (2),  (3),  and  (4)  are 


obtained, 


V  «  -U12  fps,  V  -  0,  and  V  «  -212  fps. 

|  f  | 


From  equations  (IOC),  (15A),  and  (15B)  are  obtained, 

H  *  -O.678,  p  «  1.15  and  0  •=  -0.92  . 

1  a 

Therefore  S(B)  has  the  form, 

,a 

0/„N  *B  *B  yB 

S(B)  «  — y  +  - ; - j  ^ 

\  (0.9U  R^)  (1.1  R^ 

The  translation  equations  (16A),  (16B),  and  (19A)  give, 


-  [xg  -  3.1|l  *gl  /  I  Qj  1  t 

Tg  -  txg  +  0.283  I  , 


*B  *  yB  • 

Ualng  these  axes  the  lethal  volume  was  sketched  in  figure  5. 

In  End  Game  III  the  member  is  Included  at  90"  to  the  axis  of 
the  missile.  The  velocity  vectors  at  both  the  member  and  the 
misalle  are  aligned  with  their  respective  axes.  In  this  case, 


7  •  200  fps,  V  -  0,  and  V  *  -300  fps 
x  1  » 


LETHAL  volume  for  end  GAME  I 


I 


o 


u. 

5  £ 

.  </> 
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In  this  case 
the  form, 

S(B)  - 


H  -  -l.?5,  - 

»B  ,  *B _ 

^  (.776  RjP* 


1,66,  and 

ya 

"  ( .6U*  Rj)* 


2, 111,  Hence,  S(B)  has 


0 


The  translation  equations  for  this  end  game  are. 


4  ‘  t*B  -  2*S8  *B]  / ,  ",  I  , 

y'  -  [Xg  *  O.U.  *B]  /  |  ^  | 
and  -  yB  . 


Using  these  axes  the  lethal  volume  was  sketched  in  figure  6. 

Notice  the  marked  difference  in  the  characteristics  of  the 
lethal  volume  in  End  Game  III  as  compared  to  either  End  Game  I  or 
End  Game  II.  Figure  7  shows  the  transition  in  five  stages.  The 
figures  show  End  Game  I  for  various  values  of  cos  0. 

First  consider  the  cylinder  formed  by  the  lines  Joining  the 
bases  of  the  two  cones  in  Case  1.  As  the  angle  0  increases,  the 
apex  of  the  cone  contained  in  the  cylinder  (Case  2)  approaches  the 
aide  of  the  cylinder  until  the  angle  9  equals  the  semi -cone  angle  a 
(Case  3)  .  As  9  becomes  greater  than  or  the  cone  breaks  through  the 
aide  of  the  cylinder  (Case  4)  .  At  0  *  90*  the  cone  becomes  complete 
ly  outside  the  cylinder.  (Case  5)  . 

As  0  becomes  still  larger  the  cone  originally  outside  the 
cylinder  starts  to  recede  into  the  cylinder.  At  0  ■=  180“  the 
volume  is  exactly  as  it  appears  in  Case  1,  except  the  cone  outside 
the  cylinder  at  9  «  0°  is  now  inside  the  cylinder  and  vice-versa. 

The  condition  for  a  cone  to  be  completely  contained  in  the 
cylinder  when  V  and  V  are  both  zero,  is 

7  » 


or  >  fl 
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BC  .Ml.  I  >  — !.J.l 

I  Bel  *  ./vTT'v* 

K  s 

is  required  condition. 

Following  a  similar  line  of  reasoning  gives  the 

when  either  V  or  V  is  not  zero  as 

r  * 


c 


(20A) 


(20B) 


condition 


Is  I  ,rol 


n — 7 

/  c  +  d 


(20C) 


CHAPTER  III 


MEASURES  OF  EFFECTIVENESS 

A.  INTRODUCTION 

In  this  section  expressions  for  the  two  measures  of 
effectiveness  discussed  in  the  introduction,  the  size  and 
probability  content  of  the  level  volume,  will  be  developed. 

The  scope  of  this  chapter  will  be  limited  by  assuming  a  constant 
rod  velocity  and  by  considering  the  target  to  be  a  single  linear 
member.  In  addition  the  results  of  this  chapier  are  applicable 
only  when  equations  (203)  or  (20C)  hold;  i.e.,  only  when  one  of 
the  cones  is  completely  contained  in  the  cylinder.  The  latter 
restriction  is  not  a  serious  limitation  since  the  vast  majority 
of  end-games  satisfy  these  inequalities. 

B.  THE  SIZE  OF  THE  LETHAL  VOLUME 

From  a  computational  standpoint,  the  size  of  the  lethal 
volume  is  a  very  efficient  measure  of  effectiveness .  In  this 
section,  expressions  will  be  given  with  which  the  size  of  the 

lethal  volume  can  be  evaluated  when  both  V  and  V  are  zero  and 

?  * 

when  either  V  or  V  is  not  zero. 

r  * 

Consider  the  lethal  volume  in  figure  4.  The  lethal  volume 
is  composed  of  segments  of  length  !  BCj  joining  the  surfaces  of 
the  two  cones .  The  volume  contained  in  one  of  the  cones  is  seen 
to  be  interior  to  the  lethal  volume  while  the  volume  contained  in 
the  other  cone  is  exterior  to  the  lethal  volume.  Since  these  cones 
are  identical  in  size,  the  volume  contained  in  each  is  the  same. 

The  size  of  the  lethal  volume  can,  therefore,  be  gotten  by  finding 
the  volume  of  the  skewed  elliptic  cylinder  obtained  by  joining  the 
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bases  of  the  two  cones.  A  similar  condition  is  seen  to  exist 
whenever  one  of  the  cones  is  completely  contained  in  the  cylinder; 
i,e«,  whenever  equation  (20b)  or  (20C)  holds. 

The  volume  of  the  skewed  elliptic  cylinder  is  given  by, 

V  -  tt  ab  1 1  1 1  e  1 

where, 

a  and  b  are  the  semimajor  and  semiminor  axes  of  the 
ellipse,  L  is  the  vector  connecting  the  center  of  the 
ellipses  at  the  ends  of  the  cylinder,  c  is  the  cosine 
of  the  angle  between  L  and  the  vector  perpendicular  to 
the  surface  of  the  ellipse.  This  fact  will  now  be 
related  to  the  size  of  the  lethal  volume. 

Let  V  denote  the  lethal  volume  and  |  Vl  denote  the  size  of  this 
lethal  volume.  When  Vy  *  Vz  *  0,  BC  corresponds  to  the  vector  L, 

j l  corresponds  to  e,  and  a  **  b  «  Rl*  Hence, 

l 

I  v-}  •  tt  |  BC  | !  i 
subject  to,  y 

i  >_  LlI 

1  +  V* 

When  Vv  and  V_  are  not  b*oth  zero,  the  cones  revolve  about  the 
yaxis.  The  vector  ,  (15E),  is,  therefore,  perpendicular  to 
the  surface  of  the  ellipse.  In  this  case,  BC  »  L,  a  ■  Rl>  and 


b  -  Rl/  ;  hence, 

1TSJI1 5^1 


Subject  to, 


|  qjj’  I !  B5  |  *  d3 


(22) 


(20C) 


As  an  example,  consider  one  of  the  lethal  volumes  of  the 


preceding  chapter.  For  End  Game  II,  equation  (20C)  holds,  and, 

ttr!  I  BC  1  a  . 

|  V  !  -  — -  (0.872)  -  0.825  ttRI  1  EC"  |  , 

since^ 

'  I  BC  |  [^  ,  0,  [.965,  0,  .266]  -  0.e72  1  53  I  . 

IS! 

An  important  property  of  the  size  of  the  lethal  volume  can 
now  be  stated: 

For  all  end  games  where  the  inequality  (20B)  or  the 
inequality  (2QC)  holds,  the  size  of  the  lethal  volume 
for  cuts  of  a  single  linear  structural  member  is  a  maximum 


when  the  velocity  vectors  of  the  target  and  the  missile 

are  parallel,  (V  ■  V  -  0  ),  and  the  axis  of  both  the 

r  * 

target  and  the  missile  lie  along  these  vectors  I  •  1). 

This  follows  directly  from  equations  (21)  and  (22)  since 
and  the  maximum  value  of  1  &  I  **1. 

THE  PROBABILITY  CONTENT  OF  THE  LETHAL  VOLUME 


1.  Introduction 

An  approach  similar  to  that  used  to  calculate  the  volume  1  V  | 
can  also  be  used  to  find  the  probability  content  of  V.  That  is, 
the  probability  content  of  V  can  be  considered  to  be  equal  to  the 
probability  content  of  an  elliptic  cylinder  plus  the  probability 
content  of  an  elliptic  cone  less  the  probability  content  of  an 
identical  cone.  Unfortunately,  the  probability  content  of  the  two 
cones  are  not  usually  equal;  hence,  they  do  not  cancel  as  they  did 


33 

in  the  calculation  of  V  .  In  general,  the  probability  content 
of  V,  ,  is  given  by, 

P  *=  P  +  P  -  P 

.  k  1  9  3 

where, 

P  is  the  probability  content  of  the  elliptic  cylinder, 

P  is  the  probability  content  of  the  cone  interior  to  the 
lethal  volume, 

and  P  is  the  probability  content  of  the  cone  exterior  to  the 

9 

lethal  volume. 

2 .  Background 

Figure  8  shows  a  skewed  elliptic  cylinder  and  a  right  elliptic 

cone  together  with  a  coordinate  system  (t#,  u#,  v7).  The  coordinate 

system  is  chosen  such  that  the  t/axis  lies  along  the  axis  of  the 

cone.  The  u#and  v^xes  are  chosen  to  be  parallel  to  the  major  and 

minor  axes,  respectively,  of  the  elliptic  base  of  the  cone.  In 

addition,  the  base  of  the  cylinder  is  assumed  to  be  parallel  to 

the  u^v7  plane.  The  ba^e  of  the  cylinder  nearest  the  origin  is 

assumed  to  be  centered  at  the  point  t,  '  and  the  apex  of  the  cone 

o 

nearest  the  origin  is  assumed  tc  be  located  at  t  7  .  Also  shown 
is  the  mean  (t#  ,u7  ,  v#)  of  an  arbitrary  frequency  function 
f(t#,  u*,  v7). 

Expressions  for  the  probability  content  of  these  two 
geometrical  figures  under  the  frequency  function  f(t#,  ,  v7) 

will  now  be  developed.  In  general  the  probability  content  of  a 
volume  under  the  frequency  function  f(t  ,  u  ,  ▼  )  is  given  by, 

P  -  JJJ*  f(t7,  u7,  v#)  dt7du7dv7  . 


Vol 
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The  limits  of  integration  are  easier  to  visualize  if  the 

origin  is  transformed  to  the  mean  (£*,  u*,  v  )  ,  so  let 

t  -  t'  -  t'  , 

»  -/ 

U  «  U  -  U  , 

and  v  *  v*  -  v‘  • 

The  expression  for  P  then  becomes 


p  -  Jjj  f(t,  u,  v)  dtdudv  , 

Vol 

since  the  Jacobian  of  the  transformation  is  one. 


First  consider  the  cylinder.  The  axis  of  the  cylinder  is 
formed  by  the  vector 

I-!TI[e  ,  *  ,  c  ] 

19  9  _ 

where,®  ,  *  ,  and  *  are  the  direction  cosines  of  L  with  respect 

i  3  3 

.  i  t  /. 

to  the  v*  f  u  p  v  )  coordinate  system.  The  semi -minor  and  semi- 
major  axes  of  the  elliptic  cross  sections  of  the  cylinder  are 
denoted  by  a  and  b  respectively.  The  equation  of  this  cross 
section  is  given  by, 


[\u4c(t)3a  +  [v-h(t)]3  ^ 


where. 


k(t)  and  h(t)  are  the  u  and  v  coordinates,  respectively, 
o€  the  center  of  the  cross  section  taken  at  the  point  t. 


By  similar  triangles  it  is  seen  that 
k(o)  ■  r-  ^ '  -  t  ')  -  u7 


and  that,  l  0 

k( t)  -  t  +  k(o)  -  -*■  (t  ♦  '  t  #)  *  u' 

CEO 
1  1 

Analogously,  *  .  . 

h(t)  -  ^  (t  ♦  t'  -  t  )  -  v' 
l 

Solving  equation  (23)  for  v  gives. 


v(t,  u)  -  h(t)  ±  -  7aS  -  Cu-k(t)] 

3 
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The  probability  content  of  the  cylinder  then  can  be  expressed 


as. 


where, 


and 


P 


f(t,  u,  v)  dtdudv 


u  -in*  4 1 '  - t7  , 

L  -  t  '  -  £7  , 

n  o  * 

U  -  k(t)  +  a  . 

L13  "  “  a  * 

U  -  h(t)  +  -  Vaa  -  [u-k(t)]a  , 

L  -  h(t)  -  ^  -  Cu-k(t)]3  « 

13  tt 


(2U) 


A  similar  approach  can  be  used  in  deriving  the  probability 
content  of  the  cone.  In  this  case, 

k(t)  -  k(tQ  -  tx  )  -  -vl ' 

h(t)  »  h(tQ  -  )  -  -v7 

and  are,,  hence,  constants.  However,  a(t)  and  b(t),  the  semi -major 
and  semi -minor  axes  of  the  elliptic  cross  sections  of  the  cones, 
vary  with  t.  If  the  height  of  the  cone  is  designated  by  y  ,  then 


or  »(t)  -  £  |  (t  ♦  t'  -  ts)  |  • 

b(t)  -S|  (t'-t/)!. 


Analogously, 


37 


For  the  cone,  then 

v(t,  u)  *  -V/  ±  ^(t)  -  [u+’J/]*' 

The  probability  content  of  the  cone,  P  ,  then  can  be  expressed 


where. 


u 

u 

u 

r 81 

r  22 

r  as 

p  -  J 

J 

i 

i 

8  L 

JL 

*■  L 

21 

22 

23 

dtdudv 


/  _  • 

u  ■  v  + 1  +  t  . 

*11 


L  «  t'  -  t' 
21 


+  a(t) 


L  -  -u  -  a(t ) 
S8 


U  «  -v*  +  ^7t4  «/a2(t)  -  [u  +  u'r 


“id  L  «  -v*  -  -  Cu  +  u'f  * 

*a  a(tj 

3.  An  Expression  for  the  Probability  Content  of  the  Lethal  Volume 

The  results  of  the  latter  derivations,  together  with  the 
results  of  Chapter  II,  will  now  be  used  to  derive  an  expression 
for  the  probability  content  of  the  lethal  volume. 

In  Chapter  II,  the  feasibility  of  generating  an  orthogonal 

*  •  t  • 

coordinate  system  Xg,  y^,  Zg  with  the  following  properties  was 
demonstrated: 

(a)  The  Xq  axis  was  shown  to  be  parallel  to  the  minor 
axis  of  the  cones  composing  the  lethal  volume. 

(b)  The  y*  axis  was  shown  to  be  parallel  to  the  axis  of 
this  cone  and  is  perpendicular  to  the  base  of  the  cone. 

(c)  The  Zg  axis  was  shown  to  be  parallel  to  the  major  axis 
of  this  ellipse. 
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A  lethal  volume  is  shown  in  figure  9,  together  with  such  a 
coordinate  system  constructed  with  its  origin  at  the  apex  of  the 
interior  cone.  This  coordinate  system  meets  all  the  requirements 
specified  for  the  (t#f  U#,  v *)  coordinate  system. 

Recal.'  that  the  height  of  the  cone  was  shown  to  be  -g-  ,  that 
the  major  axis  of  the  base  of  the  cone  was  shown  to  be  Rg,  and  that 

h 

the  minor  axis  of  the  base  of  the  cone  was  shown  to  be  — 

Furthermore,  if  the  vector  BC  is  restricted  to  be  directed  in 
the  positive  y^  direction  then  Lhe  vector  PC  corresponds  to  the 
vector  L.  This  restriction  does  not  limit  the  applicability  of  the 
derivation  in  any  way  since  the  choice  of  points  B  and  C  and  hence 
direction  of  the  vector  BC  is  completely  arbitrary. 

In  Chapter  III,  EC  =  1  BC|  (l  t  l  ) ,  was  given  relative  to 

l'  s'  a 

the  (x,  y,  z )  coordinate  system.  To  be  applicable  here,  BC  must  be 
transformed  to  the  (Xg,  y^,  Zg  )  system  by  equations  (16A) ,  (16B), 
and  (16C) .  The  transformed  and  restricted  vector  BC  will  be  denoted 
by  BC  *  |  BC  {  t  l'  )•  Then  «  f  c  f  e  corresponds  to  t*  t‘ f 

and  -^respectively . 

Also  shown  in  figure  9,  is  the  location  of  the  mean  y^, 

*g  )  of  the  frequency  function  f(Xg#,  y^*,  Zg#)  .  As  before  the 
origin  is  translated  to  the  mean  by  letting, 

t  -  yB'  -  7j'  , 


*  1  *B  "  *B 


and  ▼  • 


*B  "  *B 


lSce  equation  (17),  u.  20. 


y  * 


Figure  9.  COORDINATES  FOR  DETERMINING  THE  PROBABILITY  CONTE'  f  OETHELLTHAl  VOLUME 
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The  expressions  for  the  probability  content,  Y^,  of  the 
cylinder  joining  the  bases  of  the  two  cones  can  now  be  easily 
obtained.  Note  that  tQ,  the  nearest  point  to  the  origin  of  the 

h 

cylinder  is  given  by  the  height  of  the  cone  -j-.  Using  equation 


(22) ,  Pj ,  is  given  by, 

U  U  U 
11  12  13 


f(t,  u,  v)  dtdudv 


where, 


K  K  K 

li  la  13 


hi  *  t  ‘  yB  » 

*  *i  ♦  k(t)  . 

\a  “  *4.  "  * 

U  •  h(t)  +  i  M  , 


Lis  *  “  c  , 

1#  R, 

k(t)  "  ^  (t  +  yB'  -  j)  -  zg'  ,  (26a) 

« 

i#  Rt 

and  h(t)  *  p  (t  +  yB'  -  -J)  -  Xg'  .  (26b) 

a 

The  probability  content  of  the  cones  can  be  found  in  a 
similar  fashion.  For  the  interior  cone,  the  distance  from  the 
origin  of  the  (  xR/ ,  y^' ,  7.^  )  system  to  the  apex  of  the  cone, 
t  ,  is  zero.  The  probability  content  of  the  interior  cone,  P2, 


is  given  by. 


u  u  u 

ai  as  as 


P#  »  J  J  J  f(t,  u,  v)  dtdudv 


L  L  L 

81  32  83 
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where, 


u  - 

2  i 

*L 

_ t 

T  " 

yB 

L  ■ 

21 

t 

u  - 

32 

TB 

+  a(t)  , 

L  - 

22 

V 

-  a(t)  1 

U  - 

29 

(t)  -  Cu4Zf.#]2 

L  • 
as 

b(t)  /a 
a(t) 

(t)  -  [u+z^'f 

a(t)  ■ 

■  d  1 

»b'  +  t] 

1 

V't)  > 

■S> 

H  *  *1 

• 

The  form  developed  for  the  probability  content  of  the 
interior  cone  is  sufficiently  general  to  be  used  to  express  the 
probability  content  of  the  exterior  cone.  For  the'  exterior  cone, 
the  distsnce  from  the  origin  of  the  )  system  to 

the  apex  of  the  cone  is  given  by, 

t,  -  !  BO  M, 

The  quantities  k(t)  and  h(t)  are  again  constants  with, 

k(t)  -  k(  IBS |  1/  -5^')  -IWU/  , 

and  h(t)  -  h(  1 1!  I  '  -  x^)  •  !  55  1  ^  *  -  xj*  . 

The  probability  content  of  the  exterior  cone,  P3,  is  given 

b>.  nun 


?a  "J 


U  U 

_  32 


L  I* 

Si  S3 


f  3 

J  f(t,  u,  v) 


dtdudv 


US1  "  (TT5  *  I  K  I  '  -  yB' 

Uaa  “  i  *  *s  '  *  *g*  *  a(t) 

K.  -iBOM  '  -iT'  -  aft) 


and  b(t)  ■ 


!  BG  |  ^  7  -  Xg#  +  d  Va2  (t)  -  [u+Zg*  I  HJ  I  ^  'T  f 

I  BC  •  ^  -  Xg  -  d  -/aa(t)  -  [u+z^*  !  BG  I  7]3  , 

d  I  t  -(^  ♦  I  K|  -  yjj')  |  , 

1 1 1  -(-^  ♦  1  BC  I  i  '  -  )  ! 


After  evaluating  the  integrals  for  P^,  P2>  and  P3,  the 

probability  content  of  lethal  volume  P^,  is  given  by, 

P  -  P  +  P  -  P 

t  i  t  » 

The  probability  content  formulas  derived  thus  far  are  applicable 
only  when  Vy  or  Vz  is  not  zero  and  when,  in  addition,  the  inequality 
(20C)  holds. 

The  equations  for  the  probability  content  of  the  lethal  volume 
when  both  Vy  and  V2  are  zero  follow  in  a  analogous  fashion.  In  this 
caae,  the  coordinate  system  (x,  y,  z)  is  equivalent  to  the  (t#,  u#, 

i  * 

v  )  system.  If  the  mean  of  the  frequency  function  f(x,  y,  z)  ia 
(x,  y,  z),  the  substitution, 

t  *  X  -  X  , 

u  -  y  -  y  , 

and  v  -  z  -  z  , 

will  translate  the  origin  to  the  mean. 

2 

Recall  that  in  thia  case  the  cones  are  right  circular  cones  with 


?See  equation  (13),  p.  18 
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V 

heights  |  J  .  Since  the  cones  arc  circular,  a(t)  equals  b(t) 
and  in  the  cylinder  n  -  b  =  Rj  . 

Again  the  vector  BC  must  be  restricted  before  it  will 
correspond  to  the  vector  L.  If  BC  is  oriented  in  the  direction 
of  positive  x  then  BC  ■  t  BC 1  {  l  .  I  ,  t  )  corresponds  to  the 

1  *  2  ’  3 

vector  L.  BC  is  relative  to  the  (x,  y,  z)  coordinate  system  and 
&  t  l  and  ^corresponds  to  e  e  and  e  respectively.  Equations 

18  3  1  '  3  3 

(2A)  and  (25)  then  can  be  used  for  the  probability  content  of  the 
lethal  volume  in  this  case  with, 

L  «  BC 

»  V  *3  “V  V  respectively 
tf ,  u',  v/  *■  x,  y,  a  respectively 

a  «  b  «= 

‘  ■  0  (interior  cone) 
t  #  »  !  I  i  (exterior  cone) 


4.  Evaluation  of  the  Integrals 

A  scheme  for  both  hand  ev.  xuation  and  computer  evaluation 
of  the  integrals  Pi,  P2,  and  P3  will  now  be  discussed.  The  scope 
of  this  section  will  be  limited  to  the  evaluation  of  the  integrals 


with  only  the  trivariate  normal  frequency  function, 

.  3 


r(t,u,T)  - 


-Cp  *p 


a  o  a 

t  «  V 


The  integral  for  the  probability  content  of  the  cylinder  or 


the  cones  over  the  trivariate  normal  distribution  f(t,  u,  v) , 


given  above,  can  be  expressed  as, 
U  U  U 

p  "  I  J  J'  f(t»  u,  v)  dt 

v  T  T  w  * 


dt  dudv 


L  L  L 
las 


where, 


U  and  L  are  constants 

l  i 

-  k(t)  +  a(t)  , 


L#  -  k(t)  -  a(t)  , 

U#  -  h(t)  +  ||||  J ?[t)  -  [u-k(t)]“ 
and  Ls  -  h(t)  -  ^a2(t)  -  [u-k(t)]‘ 


To  evaluate  this  integral,  first  make  the  normalizing 


transformation, 


t  «=  t/a 

'  t 

u'  »  n/o 

« 

▼  #  ■  v/ct 


Equation  (27)  becomes, 

.  U  '  u  •  U  ' 


p  *  <55f^a  J.  I.  .  J.  r(t<>  u'>  T'> 


dtdudv 


where, 


L  •  L  /  L  ' 
las 


U  '  «  U  /a  . 

i  i  t  ' 

l»  '  "  \/UK  » 

O.'-kV)  ♦  a'(t')  , 


L  '  -  k'(t')  -  a'(t')  , 


0.  '  -  fc'(t')  *  SJli  ,/  [a'(t')]1  -  [u'-k'(t')]a  , 

s  a  (t  ) 

L  '  -  h'(t. ')  -  ,/  [a'(t')]a  -  [u'-k'(t')]s  , 

k*  »  k/o>  , 

h  ■  h/a  # 

a  -  a/a  , 

/  v /» 

and  b  -  b  /a  , 
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Since  this  transformation  maps  ellipses  into  ellipses  and 
straight  lines  into  straight  lines,  the  problem  is  still  to 
evaluate  the  integral  Y  over  a  right  elliptic  cone  and  skewed 
elliptic  cylinder. 

A  computer  program  suitable  for  the  evaluation  of  P  is 

3 

readily  available  for  the  IBM  7090  computer.  This  routine 
computes  the  integral,  P,  using  a  Hasting's  approximation^  for 
the  inner  integral  then  sums  by  n  step  by  step  process. 

The  integral  P,  can  be  evaluated  in  terms  of  tabulated 
functions  only  if  k #(t /),h#(t#),  are  constants  and  a*  (t)  and 
b*(t)  are  equal.  Since  a^t/)  and  b  (t  )are  always  different  for 
the  cones  when  Vy  or  Vz  is  not  zero,  numerical  integration  must 
usually  be  used  to  evaluate  the  probability  content  of  the  cones. 
With  the  cylinder,  a(t)  and  b(t)  are  always  constants  and  when 
either  Vy  or  V2  is  not  zero,  equations  (26A)  and  (26B)  show  that 
k*(t#)and  h#(t/)  are  also  constants  if  *  •  0  «  The 

corresponding  condition  for  the  case  when  both  V^  and  Vz  are  zero  is 

h  »  X  ■  0*  att  t 

When  a*(t#),  b#(t*),  k#(t#)  and  b  are  constants  then 


^Gittleman,  Ruth,  Floating  Point  (n)  Variate  probability 
Integral,  IBM  SHARE  ROUTINE  1384,  November  2,  1962. 


^Hastings,  Cecil,  Jr.,  Approximation  for  Digital  Computers, 
Princeton  University  Press.  1955,  p.  167. 


r 
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the  integral  >  can  be  written  as 


/  >  3 


Of 


.-W 


k  +a  nh  +  ~  /a  -  (u  -k  ) 
a 

i  -i(u#+v/)2du#dv# 
/  i  v  e  ^ 


I  - - 

k  -a  h - -  .'a  -  (u  >k  ) 


or 


where, 


P  ■  AO  E(a#,  b#,  k*,  h#) 

An  -  n(u  ')  -  n(L  ')  , 

°jio>  •  r e4t* dt/ 


(28A) 

(28  b) 

(28C) 


and 


a 


Q(X)  is  the  integral  of  normal  distribution  function  and  is 
tabulated  in  numerous  books. 

E(a',  b',  kJ  \  h#)  is  the  integral  of  the  bivariate  normal 
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distribution  over  an  offset  ellipse.  This  function  is  tabulated 
in  at  least  three  places.^’  ^  The  notation  was  chosen  to 
conform  to  that  of  Mr.  Builte.8  J.  H.  Cadwcli,^  has  shown  that 


E(a  ,  b  ,  k  ,  h  )  «  --  - 


2a  V 


4 


re.*  “■)] 


exp 


(20) 


a 


'  ♦!»  t# 


Mr.  Cadwell  showed  that  the  error  from  using  tiiis  approximation 
is  less  than  .002  if  a  is  less  than  0.5.  He  also  gives 
approximations  which  can  be  used  for  a#less  than  1.5  and  for  a 
less  than  2.5.  However,  with  comparable  accuracy  these  latter 
approximations  are  rather  cumbersome. 


5Rodden,  J.  J.  and  G.  W.  Rosenthal,  Tables  of  the  Integral  of 
the  Elliptical  Bivariate  Normal  Distribution  Over  Offset  Circles, 
Lockheed  LMSD  -  900619,  May  1,  1961 . 


^DiDonato,  A.  R.  and  M.  P.  Jarnagin,  Integration  of  the 
General  Bivariate  Gaussian  Distribution  Over  an  Offset  Ellipse. 
Naval  Weapons  Laboratory  Report  NWL  1710,  August  11,  1960. 


^Builte,  A.  G.,  A  Building  Block  Technique  for  the  Statement 
and  Solutions  of  the  Problems  Involving  the Bivariate  Normal 
Distribution.  Master's  Thesis,  School  of  Engineering,  Air  University, 
August  1962,  p.  119,  pp.  54-50. 


8Ibid. 


^Cadwell,  J.  H.,  "An  Approximation  to  the  Integral  of  the 
Circular  Gaussian  Distribution  Over  an  Offset  Ellipse",  Mathematics 
of  Computation.  January  1964,  Volume  18,  No.  85,  pp.  106-113. 


In  the  special  case  where  a'  *  b#,  the  circular  coverage 
function,  p(a#,  r) ,  of  H.  H.  Germond  is  applicable.*®  In  this 


case. 


E(o',  b',  k',  b)  -  E(a',  b,  k',  h')  -  p(a',  r) 

where ,  ,»  J>  * 

r  ■  h  ♦  k  • 

An  extensive  tabulation  of  the  o  function  where, 

Q(a',  r)  -  1  -  p(a',  r) 

exists  and  offers  the  best  source  of  the  values  of  the  circular 
coverage  function. ** 

Roger  Snow  showed  that  the  circular  coverage  function, 
p(a#,r  )  can  also  be  used  when  a  /  b  if  k  ■  h  In  this 


•  * 


case, 


E(a',  b',  0,  0)  -  q(a',  b#)  -  p(^-^-A^-)  -  p[  ' *  g*  ■  ,  &  g*"'']* 
If  both  a7®  b*  and  k#  •  h*  ■  0,  then  E(a#,  a#,  0,  0)  *  l-«^a 


/  / 


i.e,,  E(a  ,  a  ,  0,0)  is  the  integral  of  the  well-known  Rayleigh 
distribution.  In  this  case  the  tables  of  e*"*  can  be  used  to 
evaluate  the  integral.  If  the  accuracy  requirements  are  not 


l0Germond,  H.  H.,  The  Circular,  i 
Corporation  RM  330,  January  26,  1950 


l»  Rand 


**Marcum,  J.  I.,  Tables  of  Q  Functions,  Rand  Corporation 
RM  339,  January  1950. 


l^Snow,  Roger,  Some  Characteristics  of  the  Elliptic  Gaussian 
Distribution.  Rand  Corpi ration  RM  2765-FR,  September  1961. 


^National  Bureau  of  Standards,  Tables  of  the  Exponential 
Function  ex,  U.  S.  Government  Printing  Office.  1947. 


to  be  stringent,  special  graph  paper  can  be  used  to  evaluate 


E(a',  a',  0,  0>.14 

If  k  (t')andh  (tj  or  a#(t*)  and  b#(t/)are  not  constants  the 
integral  can  be  evaluated  by  the  following  procedure: 

(1)  For  a  given  t  #«  *  evaluate  k/(tj,)and  h#(tj  ') 

or  a^t^  /)  and  b/(tj  /). 

(2)  Find  E  ta#(t^  b^t^  #),  Y.' t  h*  ]  from  one  of  the 
references  given  or  from  equation  (29) . 

(3)  From  the  normal  probability  tables  find 


obabiiity  tables  tir 
v  .  1  . 73  t  r 
)  *  (v=y  e  5  I 


(A)  From  the  product, 

^  ^  ^  ) 

This  determines  a  point  of  the  frequency  function, p(t  ). 

(5)  Repeating  this  procedure,  with  t  Varying  between  the 

lower  limit  L  and  upper  limit  U  ,  determines  the 
1  1 
entire  frequency  function. 

(6)  Numerical  integration  of  this  frequency  function  yields 
the  probability  content  P. 

5.  Example 

As  an  example  of  the  above  procedure  the  probability  content 
of  the  lethal  volume,  will  now  be  computed  for  a  special  case  of 
End  Game  I.  In  End  Game  I  if  f  *1,  and  t  ^  *0,  then  an 

1  S3 

antiparallel  intercept  exists.  In  this  case  a  ■  b  ■  and  the 
cylinder  is  right  circular.  If  the  mean  of  the  distribution  is 
assumed  to  be  located  at  the  base  of  the  interior  cone,  (see  figure 
10);  then,  k#»  h*0»  and  equations  (28A)  and  (30)  hold. 


l^Burkc,  T.  Finley,  New  Graph  Paper  for  Circular  Normal 


Distributions ,  Rand  Corporation  RM  3292 -PR,  September  1962. 
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To  fix  the  dimensions  of  Che  lechal  volume,  lec  ’  BC 1  «*  R,  and 
V 

«r  ■  To  iix  Che  discribuCion,  let  3  ■  cr  »  a  *  ft.  # 

VR  *  r  *  V3 

Tor  Che  cylinder,  Chen  U  -  3.0  and  L  *  0^0.  EquaCions 

11  11 

(28B)  and  (?.8C)  Chen  yield, 

l  t* 

An  .  — i-r  I  e"  T  dt  -  O.U987 

<2">iJ0 


Hence 


From  equacion  (30)  is  obcained 

E  »  1  -  -  0.98889 

P  -  AQ  E  *  0.4940 

Nexc  consider  Che  problem  of  finding 


p3* 


From  figure  10 


Che  cone  excerior  co  Che  volume  is  seen  co  be  orienced  wich  ics 

apex  locacion  ac  t#  *  +1.5  ■  L  .  Since  V  /VD  ■  i,  U  ■  L 

*  ai  *  «  ai  ai 

♦  (i)  Rj/ot  -  3.0. 

Following  Che  procedure  ouc lined  in  Che  lasc  sec  cion,  Table  2 
can  now  be  generaced. 


Table  2.  Calculation  of  : 

P  a 

JL 

•  V)  1  - 

8  1  -V 

*,<*') 

1.5 

0  0 

.12952 

0 

1.6 

0.2  .01980 

.11092 

.00220 

1.7 

• 

0.4  .07688 

«  • 

.09405 

• 

.00723 

• 

• 

• 

3.0 

«  « 

•  • 

3.0  .98887 

• 

• 

.00443 

• 

< 

.00438 

Where 

P.(t')  -[l  -  °'iC*V)rl(s; 

? 

)  . 

normalized  lethal  volume  for  antiparallel  intercept 
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Next  p  (Ois  plotted  in  figure  11.  The  area  uader  this 

9 


curve  Is  estimated  using  one  of  numerous  techniques  available 

3.0 

P  *  f  p  (t/)  dt  ■  0.0188  . 

a  J1.5  3 


yielding. 


Pj  can  be  found  by  a  similar  procedure.  From  figure  10  the 
cone  interior  to  the  volume  is  seen  to  have  its  base  nearest  the 
mean  of  the  distribution.  For  the  interior  cone  the  function 
•V)  must  be  a  maximum,  (3),  at  t/«0  md  a  minimum,  (0), 
at  The  function  p  (t/jis  plotted  in  figure  12.  Hen 

•  3 

«  f°  pjt')  dt  -  0.285 

of  the  volume  is  given  by 

P  -  0.7602 

s 


•  J  » 

-1.5 

Now  the  probability  content 


-  P  ♦ 
l 


P  - 
* 
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CHAPTER  4 

PARAMETER  VARIATION  STUDIES 

Using  the  results  of  Chapter  3,  the  variation  in  the  size 
and  the  probability  content  of  the  lethal  volume  vith  the  model 
parameters  can  be  investigated.  In  a  typical  investigation  all 
except  one  of  the  model  inputs  are  held  constant  and  the  variation 
in  the  measures  of  effectiveness  is  expressed  as  a  function  of  the 
remaining  parameters.  The  values  assigned  to  the  inputs  held 
constant  are  then  changed  «*nd  the  variation  in  the  measures  of 
effectiveness  is  again  expressed  as  a  function  the  remaining 
parameter.  This  process  is  repeated  until  the  change  in  the 
measure  of  effectiveness  can  be  reliably  predicted  for  any  set 
of  values  assigned  to  the  Inputs  held  constant. 

Another  way  In  which  the  variation  in  the  measures  of 
effectiveness  is  sometimes  investigated  is  the  tradeoff  study. 

For  the  tradeoff  study  all  except  two  of  the  inputs  are  held 
constant.  Combinations  of  these  two  inputs  are  then  found  such 
that  the  measures  of  effectiveness  remain  constant.  Here,  a 
popular  example  is  a  set  of  iso*-klll  probability  curves.  Having 
such  curves  the  analyst  can,  if  he  has  the  cost  data,  select  the 
combination  of  the  two  variables  which  give  a  satisfactorily  high 
kill  probability  for  the  least  cost. 

The  construction  of  iso-kill  probability  curves  in  this 
fashion  may  be  extremely  difficult  unless  the  probability  of 
kill  function  can  be  solved  explicitly  for  the  variables  of 
Interest.  Since  the  kill  probability  function  cannot  usually  be 


w 


T- 
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solved  explicitly,  an  iterative  process  must  be  used.  This  can 
be  quite  expensive  if  high  accuracy  is  desired. 

One  way  this  problem  can  be  surmounted  is  to  compute  the 
probability  of  kill  as  a  function  of  the  ratio  of  two  of  the  model 
Inputs.  If  the  probability  of  kill  can  be  computed  in  this  fashion, 
then  the  tradeoffs  between  two  of  the  inputs  can  still  be  noted  and 
the  calculations  can  still  be  easily  performed. 

A  study  in  which  the  variation  of  kill  probability  is  expressed 
as  a  function  of  the  ratio  of  the  input  variables  will  not  be 
Illustrated.  To  obtain  the  ratios,  the  trivariate  normal  frequency 
function  f(t,  u,  v)  Is  restricted  to  the  spherical  normal  frequency 
function;  i.e., 

a  m  a  M  o  *  o  . 

t  I  T 

The  probability  content  of  the  lethal  volume  can  now  be 
computed  as  a  function  of  the  ratio  of  the  maximum  opening  radius 
of  the  rod,  RL,  to  or  as  a  function  of  the  ratio  of  the  effective 
length  of  the  member  I BC  l ,  to  o# 

Figure  13  shows  the  location  of  the  mean  of  the  frequency 
function  taken,  the  values  assigned  for  the  other  inputs,  and  the 
range  of  the  ratios  that  were  investigated. 

The  results  of  this  investigation  are  shown  in  figures  14  and 
15.  Figure  14  shows  the  size  of  the  lethal  volume  normalized  with 
respect  to  <7*as  a  function  of  the  ratio  Rj/<7.  Figure  15  ahows 
the  variation  In  the  probability  content  of  the  lethal  volume 
for  the  same  inputs. 

The  lethal  volume  model  can  be  similarly  used  to  determine 
the  vsriation  in  kill  probability  with  any  of  the  parameters  which 


ENO  GAME 


Figurt  13.  ASSUMPTIONS  FOR  FIGURES  14  AND  15 


THE  SIZE  OF  THE  NORMALIZED  LETHAL  VOLUME,  V{A)|/<r3 


Figure  -  14.  THE  VARIATION  IN  THE  SIZE  OF  THE  LETHAL  VOLUME 
WITH  THE  MODEL  INPUTS 
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«re  Inputs  to  the  model.  The  number  of  parameters  of  the  model 
can  be  increased  by  deriving  the  relationship  between  the  inputs 
and  other  end  game  variables.  For  example,  the  standard  deviation 
of  frequency  distribution  can  be  divided  into  guidance  and  timing 
errors.  Parameter  variation  studies  can  then  be  conducted  with 
either  the  guidance  or  timing  errors.  A  further  division  such 
as  thia  Increases  the  resolution  of  the  model.  Table  3  lists 
some  of  the  model  inputs  together  with  ways  these  inputs  might 
be  divided  to  increase  the  number  of  parameter  variation  studies 
possible  with  the  model. 

Table  3 

Division  of  the  Model  Inputs 

1.  Standard  deviation  of  the  frequency  distribution 

a .  Guidance  errors 

b.  Timing  errors 

c.  Target  maneuvers 

2.  Lethal  radius  of  the  rod 

a .  Warhead  weigh  t 

b.  Warhead  size  and  shape 

3.  Average  rod  velocity 

a.  Warhead  weight 

b.  Warhead  Size  and  shape 

c.  Engagement  altitude 

4.  Effective  length  of  target 

a .  Target  shape 

b.  Target  hardness 

c.  Closing  velocity 
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5. 


Velocity  vector  of  the  missile 

a.  Range  to  intercept 

b.  Angle  between  the  missile's 


axis  and  its  velocity  vector 


CHAPTER  5 


MODEL  IMPROVEMENTS 

INTRODUCTION 

In  this  final  section,  a  means  of  improving  the  target  and 
warhead  models  will  be  presented.  A  model  improvement  is  defined 
as  a  more  realistic  mathematical  description  of  the  actual 
performance  of  the  hardware  during  intercept. 

The  need  for  model  improvements  depends  on  the  particular 
variables  to  be  investigated,  the  purpose  of  the  study,  the  accuracy 
required,  and  the  reliability  of  the  inputs  to  the  model.  For 
example,  if  the  object  of  a  study  is  to  accurately  estimate  the 
kill  probability  of  a  system,  then  a  more  realistic  description 
of  the  target  and  warhesd  would  be  desirable.  On  the  other  hand, 
if  the  object  of  the  study  is  to  roughly  estimate  the  gain  in 
kill  probability  obtained  by  improving  the  guidance  accuracy,  then 
the  current  model  will  suffice. 

TARGET  MODEL  IMPROVEMENTS 

The  lethal  volume  associated  with  cuts  of  a  structural  member 
haa  been  derived.  The  lethal  volume  associated  with  cuts  of  a 
stick  aircraft  can  now  be  determined  by  taking  the  union  of  the 

♦ 

lethal  volume  associated  with  each  of  the  structural  members  of 
an  aircraft-  One  of  the  ways  that  this  may  be  accomplished  is 
given  below. 

Consider  the  aircraft  sketched  in  figure  16.  If  the  point 
where  the  wings  are  attached  to  the  fuselage  is  selected  as  the 
origin,  then  the  target  can  be  desciibed  with  four  vectors.  One 
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vector,  N,  has  a  magnitude  equal  to  the  length  from  0  to  N,  and  is 
directed  from  0  to  N.  Similarly,  vectors  from  T,  and  W2  are 
defined  as  the  vectors  from  0  to  T,  W^,  and  W2  respectively. 

The  total  lethal  volume  Vl»  associated  with  cuts  of  all 

four  vectors  is  then  given  by,  ^  ^  ^ 

1  a 

Figure  17  gives  the  cross-section  of  Vl  in  the  xy  plane  for 
ah  antiparallel  intercept  (End  Game  I  with  ^  •  3.)  .  This  approach 
ignores  the  cumulative  effects  which  might  occur  when  more  than 
one  structural  member  is  hit;  i.e.,  when  two  members  are  both 
subjected  to  a  sub-lethal  effect,  their  cumulative  effect  may  well 
be  lethal  to  the  aircraft.  The  changes  of  such  an  occurrence  are 
remote  with  a  continuous  rod  warhead. 

The  major  difficulty  that  arises  in  making  this  model  improve¬ 
ment  is  in  expressing  the  limits  for  the  integral. 


P 

k 


f(t,  u,  v)  dtuddv 


Ideally  the  limits  of  integration  could  be  expressed  as  functions 
of  the  wing  length  IVT"!  ,  the  wing  angle  m,  and  the  roll  of  the 
aircraft.  Such  a  derivation  would  be  an  extremely  difficult  but 
certainly  a  feasible  undertaking. 

To  better  estimate  the  desirability  of  undertaking  the  target 
model  improvement,  an  estimate  of  the  change  in  kill  probability 
due  to  adding  the  wings  to  the  target  is  needed.  If  the  vector 
NT  is  considered  to  represent  the  segment  l  BC?  for  the  single 
linear  structural  member  model  already  developed,  then  the 
consequences  of  adding  the  vectors  W*  and  TT to  the  model  can 


be  considered. 
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From  equation  (31)  it  is  seen  that  since 

V  VHTU  <V  V,)-VF 

the  probability  content  with  the  stick  aircraft  model  will  in 

general,  be  greater  than  that  obtained  by  using  the  single 

member  model.  In  figure  17,  the  cross  section  of  -  (V^  U  ) 

l  a 

has  been  shaded.  The  probability  content  of  -  (Vy  U  ) 

l  a 

represents  the  increase  in  probability  due  to  the  addition  of  the 
vings .  For  a  particular  end  game  the  probability  content  of 
Y-  -  (Vw  U  VW)  can  be  bounded  in  numerous  ways,  e.g.,  using 

l  a 

Gaussian  equal  cell  probability  paper.  From  this  Increase  in 
probability  the  analyst  can  Judge  the  importance  of  the  target 
model  Improvement  for  his  particular  application. 

WARHEAD  MODEL  IMPROVEMENTS 

One  of  the  principal  assumptions  made  before  the  expressions 
for  the  size  and  probability  content  of  the  lethal  volume  are 
derived  in  Section  D,  was  that  the  rod  travels  with  a  constant 
velocity.  In  the  sequel,  a  means  of  accounting  for  a  variable 
rod  velocity  and  the  effects  of  this  model  Improvement  will  be 
discussed. 


If  an  object  is  ejected  into  the  atmosphere  with  high 
initial  velocity,  its  velocity  will  decay  rapidly  as  a  function 
of  the  distance  traveled.  This  velocity  will  approximate  the 


equation, 

where. 


V(R)  ■  V  e'1® 
o 


(32) 


V(R)  is  velocity  of  the  rod  at  the  distance  R, 
V(o)  is  the  ejection  velocity  of  the  rod, 
and  K  is  the  decay  parameter. 
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The  parameter  K  Is  a  function  of  the  drag  encountered;  hence, 

K  varies  with  the  ballistics  of  the  rod  and  with  the  air  density. 
The  parameter  K,  then  will  not  oily  change  with  the  characteristics 
of  the  rod,  but  will  also  change  with  the  altitude  of  the  burst. 

The  time  required  for  the  rod  to  travel  a  distance  R  is 


given  by. 


'n  ‘  l  WJ  '  lTVToT  (•**; 


The  average  velocity  of  the  rod  to  a  distance  R,  the  parameter 
Vr  used  in  Chapter  II,  is  given  by. 


R  K  R  V(0 


e  -1 


When  V,  =  V-  c  0,  equation  (13)  gives 

V  * 


.  *  T  a 

S(B) 

R*  R* 


Where  S(B)  -  0  is  the  equation  of  the  surface  associated  with  cuts 
of  the  point  B. 

Substituting  for  VR  gives. 


s'(B) .  !s_  ♦  Is. _ i 


h*  if 


rK  V  V  R*-. 

bs-J 


Following  a  procedure  similar  to  that  unsed  in  Chapter  3, 
the  lethal  volume  can  be  formed  by  connecting  the  surface  defined 
by  S#(B)  "0  to  a  similar  surface  defined  by  S#(C)»0with 
parallel  line  segments.  Then  the  expressions  for  the  size  and 
probability  content  of  the  lethal  volume  can  be  developed  following 
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the  techniques  of  Chapter  4.  Again  the  major  difficulty  would  be 
in  expressing  the  limits  of  the  integration  to  obtain  the 
probability  content. 

Note  that  since  K  is  greater  than  zero  the  surface  defined 
S#(B)*0  contains  a  larger  volume  than  the  surface  defined  by 
S(B)  ■  0«  (See  figure  18)  .  Again  an  estimate  of  the  change  in 
kill  probability  due  to  the  model  improvement  is  desired. 

If  is  the  probability  content  of  the  surface  defined  by 

,  .  » 

S(B)  ■  0  (obtained  with  VR  ■  VR^)  and  *s  t*ie  Fr°bability  content 

of  the  surface  defined  by  S#(B)  *  0  (using  equation  (32)  for  VR) 

then  the  problem  is  to  estimate 


Let  the  axis  of  the  surface  defined  by  S*(B)  «=  0  be  divided 

into  n  intervals  and  let  j  »  1,  2,  --  n,  be  the  coordinate  along 

the  axis  at  the  end  of  the  n1*1  interval.  Let  R#  denote  the  distance 

J 

from  the  axis  to  the  surface  defined  by  S  (3)  at  the  point  J. 

Let  the  probability  content  of  the  cones  of  radius  R *  obtained 

with  VR  ■  VRj  be  P j .  Let  Rj  be  the  distance  from  the  Jth 

coordinate  to  the  cone  S(B)  »  0  and  let  Pj  be  the  probability 

content  of  this  cone* 

Note  that  P  *  m  P  “  Pr*  Then, 

Pc' -  P0  -  Urn  E  (P  '  -  P  )  03) 

n— •  j-1 

taking  a  small  n  will  yield  the  desired  approximation. 

A  similar  procedure  will  yield  P  P  *  and  P  *  corresponding 
to  the  probability  content  of  the  cones  and  the  cylinder.  The  same 
procedure  can  also  be  used  when  either  Vy  and  Vz  are  not  both 
aero.  The  analyst  can  again  judge  whether  or  not  this  model 


improvement  is  worthwhile. 

Incidentally,  equation  (33)  offers  an  alternate  means  of 
obtaining  the  probability  content  with  variable  rod  velocity. 
Using  equation  (33)  to  obtain  a  high  degree  of  accuracy  would 
be  rather  expensive,  however. 
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